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Abstract 

The self-similar growth-fragmentation equation describes the evolution of a 
medium in which particles grow and divide as time proceeds, with the growth 
and splitting of each particle depending only upon its size. The critical case 
of the equation, in which the growth and division rates balance one another, 
was considered by Doumic and Escobedo [11] in the homogeneous case where 
the rates do not depend on the particle size. Here, we study the general 
self-similar case, using a probabilistic approach based on Levy processes and 
positive self-similar Markov processes which also permits us to analyse quite 
general splitting rates. Whereas existence and uniqueness of the solution are 
rather easy to establish in the homogeneous case, the equation in the non- 
homogeneous case has some surprising features. In particular, using the fact 
that certain self-similar Markov processes can enter (0, oo) continuously from 
either 0 or oo, we exhibit unexpected spontaneous generation of mass in the 
solutions. 
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1. Introduction 

The growth-fragmentation equation is a linear differential equation intended to 
describe the evolution of a medium in which particles grow and split as time passes. 
It is frequently expressed in terms of the concentration of particles with size a: > 0 at 
time t, say u{t,x), as follows: 

pOO 

dtu(t,x)+dx{T{x)u{t,x))+B{x)u{t,x)= k{y,x)B{y)u{t,y)dy, (1.1) 

J X 

where T(a:) is the speed of growth of a particle with size x, B{x) the rate at which a 
particle of size x splits, and k{y,x) = k{x — y,x) twice the probability density that 
a particle with size x splits into two particles with size yx and (1 — y)x (the factor 2 
is due to the symmetry of the splitting events). This type of equation has a variety 
of applications in mathematical modeling, notably in biology where particles should 
be thought of as cells, and has motivated several works in the recent years; see, for 
example, [11], which also contains a summary of some recent literature. 
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We are interested here in the situation t(x) = B(x) = x°‘ for some a G K 

and k has the form k{y,x) = x~^ko{y/x); for these parameters, (1.1) possesses a useful 
self-similarity property. This is referred to as the critical case by Doumic and Escobedo 
[11], who studied in depth the situation when, additionally, a = 0. For our purposes, it 
will be more convenient to write the equation in weak form, as follows. For a; > 0 and 
y G [1/2,1), we write (y | x) for the pair {yx, (l — y)x}, which we view as the dislocation 
of a mass x into two smaller masses, and then for every function /: (0, oo) —>■ M, we 
set 

f{y I x) ■■= fiyx) + /((I - y)x). 

Consider test functions / G C/°(0, oo), that is, / is inhnitely differentiable with com¬ 
pact support. For a measure y on (0, oo), we write {y, f) ■= f{x) y{dx). By 

integrating (1.1), we obtain the equation 

dt{ytJ) = {yt,Cf), ( 1 . 2 ) 

where ytidx) = u{t, x)dx and the operator £ has the form 


Cf{x) = x°‘ cxf{: 


^)+ / 

Jli 


if{y I x) - f{x))K{dy)\ , 
[ 1 / 2 . 1 ) / 


(1.3) 


where 

K{dy) := koiy)dy = ko(l - y)dy , y G [1/2,1), 

is referred to as the dislocation measure. The advantage of this formulation is that 
we do not require absolute continuity of the solution yt or the dislocation measure 
K. More generally, one might also consider non-binary dislocation measures, but we 
refrain from doing so in this work in order to simplify the presentation. 

In this article, we study the equation (1.2) for operators of the form 


Baf{x) ■■= x°‘(ax'^f"{x) + bxf{x) + f {f{y \ x) - f{x) + xf{x){l - y))K{dy)] , 
\ “'[ 1 / 2 . 1 ) / 

(1.4) 

where a > 0, & G M, and we now only assume that the measure K satisfies the weaker 
requirement 

[ - y)‘^K{dy) < oo. (1.5) 

-'[i/a.i) 

Our notion of a solution of (1.2) is a collection of locally hnite measures {yt)t>o on 
(0,oo) such that, for every / G C/°(0,oo) and t > 0, there is the identity 

{l^tJ) = {yo,f)+ [ {ys,Cf)ds. 

Jo 

(This requires implicitly that s i-G {ys, Cf) be a well-defined, locally integrable function, 
and in particular the family {yt)t>o is then vaguely continuous). 

We offer a comparison between the original operator (1.3) and our operator (1.4). 
Besides the appearance of a second order derivative, there is a new term xf'{x)(l — y) 
in the integral in (1.4). The latter should be interpreted as an additional growth term 
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which, in some sense, balances the accumulation of small dislocations. We stress that 
(1.5) is the necessary and sufficient condition for (1.4) to be well-defined, and that 
when the measure K is finite (or at least fulfills ~ every 

operator of the form (1.3) can also be expressed in the form (1.4). Our motivation for 
considering this more general setting stems from the recent work [3], in which a new 
class of growth-fragmentation stochastic processes is constructed such that, loosely 
speaking, the strong rates of dislocation that would instantaneously shatter the entire 
mass can be somehow compensated by an intense growth; the dislocation measure 
associated with such a fragmentation process need only satisfy (1.5). 

In short, the purpose of this work is to demonstrate the usefulness of some proba¬ 
bilistic methods for the study of these critical growth-fragmentation equations. More 
precisely, we shall see that solutions to (1.2) for £ = £„ can be related to the 
one-dimensional distributions of certain self-similar Markov processes, and this will 
enable us to reveal some rather unexpected features of the former. Although in the 
homogeneous case a = 0, we establish existence and uniqueness of the solution in full 
generality, this feature is lost for a ^ 0. In particular, we shall see that under a fairly 
general assumption on the parameters of the model, the critical growth-fragmentation 
equation permits spontaneous generation, i.e. there exist non-degenerate solutions 
starting from the null initial condition. 

We need some notation before describing more precisely our main results. We first 
introduce the function k: [0, oo) —> (—oo, oo] which plays a major role in our approach 
and is given by: 


K{q) := aq^(b - a)q{I - - 1q{l - y)) K(dy), g > 0. (1.6) 

J [1/2,1) 

There are two principal ‘Malthusian hypotheses’ which we will require when a 0: 

(M+) infq>o K{q) < 0. 

(M_) There exist 0 < a;_ < w+ and e > 0 such that k{oj_) = k{oj+) = 0 and 
k{uj- — e) < oo. 

We now summarise our main results, deferring their proofs to the body of the article. 

• For a = 0, the equation (1.2) with operator (1.4) and initial condition /io = (5i 
has a unique solution. 

• For a < 0: suppose that (M+) holds. Then, the equation (1.2) with operator 
(1.4) has a solution with initial condition /io = (5i. There exists further a non¬ 
degenerate solution started from /io = 0; in particular, uniqueness fails. 

• For a > 0: if (M+) holds, then the equation (1.2) with operator (1.4) has a 
solution with initial condition pto = Si. If (M_) holds, then there also exists a 
non-degenerate solution started from /io = 0; again, uniqueness fails. 

We shall also observe that under essentially the converse assumption to (Af_|_), 
namely that infg>o K{q) > 0, the particle system that corresponds to the stochastic 
version of the model may explode in finite time almost surely. This is a strong indication 
that (1.2) should have no global solution in the latter case. 
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The rest of this article is organized as follows. The next section provides brief prelim¬ 
inaries on the function k and the use of the Mellin transform in the study of growth- 
fragmentation equations. Section 3 is devoted to the homogeneous case a = 0, and then 
the general self-similar case a 7 ^ 0 is presented in section 4. In section 5, we investigate 
a stochastic model related to the growth-fragmentation equation, to demonstrate that 
explosion may occur when the Malthusian hypothesis fails. Finally, in section 6 , we 
briefly discuss another interpretation of the growth-fragmentation equation in terms of 
branching particle systems and many-to-one formulas, placing the results of sections 3 
and 4 in context. 

2. The Mellin transform and the growth-fragmentation equation 

We observe first that, for any a G R, the operator £„ fulfills a self-similarity 
property. Specifically, for every c > 0, if we denote by (fidx) = cx the dilation function 
with factor c, then for a generic / G C“( 0 , 00 ), there is the identity 

O Vc) = C~°‘ [Caf) O (2.1) 

As a consequence, if is a solution to (1.2) for all / G C“(0, cx)) with initial 

condition and if pt denotes the image of pt by the dilation pc, then 

is a solution to (1.2) for all / G C“(0, oo) with initial condition po = Sc- For the sake of 
simplicity, we shall therefore focus on the growth-fragmentation equation with initial 
condition po = 5i, since this does not induce any loss of generality. 

Recall that the function k has been introduced in (1.6); its domain is clarified by 
the following result. 

Lemma 2.1. 

(i) For every q >0, K,{q) is well-defined with values in (— 00 , 00 ]. The function k is 
convex, and we define dom^ = {g > 0 : K{q) < 00 }. 

(a) For g > 0, K{q) < 00 if and only if /[ 1/2 i)(^ ~ yYK(dy) < 00 , and in particular 
[ 2 , 00 ) C dom K. 

(Hi) For every function f in C“(0, 00 ), £„/ is a continuous function on (0, 00 ) and 
is identically 0 in some neighborhood of 0. Furthermore, Caf{x) = o{x^~^°‘) as 
X ^ 00 for every q G dom k, and thus in particular for q = 2. 

Proof. (i“ii) First, the integral 1 ) (j/® ~ 1 + 9(1 ~ 2/)) K{dy) converges ab¬ 
solutely thanks to (1.5), since y'^ — l-\-q{l — y) = 0((1 —y)^). It follows that K{q) is 
then well-defined with values in (— 00 , 00 ) if and only if J^i /2 i)(l“ 2 /)'^-^(*l 2 /) < 00 , 
and otherwise K{q) = 00 . 

(iii) The first assertions are straightforward, and so we check only the last one. 
Take q G domK and recall from above that i)(l ~ ydK{dy) < 00 . This 
entails Ar([l/2,1 — e)) = o{e~d as e —)■ 0+. Since / has compact support in 
( 0 , 00 ), we have for x sufficiently large that Cof{x) = 1 ) fi^i^ ~ y))^i'^y)j 

and we easily conclude that Cof{x) = o{x'^). 

Doumic and Escobedo [11] studied certain growth-fragmentation equations with 
homogeneous operators given by (1.3) for a = 0, and observed that the Mellin trans¬ 
form plays an important role. In this direction, it is useful to introduce the notation 
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hq-. ( 0 , cx)) —>■ ( 0 , cx)), hq{x) = x‘^, for the power function with exponent q, and recall 
that the Mellin transform of a measurable function /: (0, oo) —> R is defined for 2 ; € C 

by 

poo 

Mf{z) := / f(x)x^-^dx 
Jo 

whenever the integral in the right-hand side converges. It follows from Lemma 2.1 that 
the Mellin transform of Caf is well defined for all z < —2 — a, or more generally for 
all z such that —z — a € (dom k)° . 

The role of k in this study stems from the following lemma, which is easily checked 
from elementary properties of the Mellin transform; see [11, §1.1] and [10, §12.3]. 

Lemma 2.2. 

(i) LetqGdomK. Then, 


(^) 7 X ^ 0. 

In particular, for a = 0, hq is an eigenfunction for Cq with eigenvalue K{q). 

(a) For every q such that q — a € (domK)° and every f € Cf°(0,oo), there is the 
identity 


M{Caf){-q) = n{q - a)Mf{-q -f a). 


3. The homogeneous case 

Throughout this section, we assume that a = 0, and refer to this case as homo¬ 
geneous. Recall that when a = 0 and the dislocation measure K fulfills the stronger 
condition 

[ -y)Kidy) < 00 , (3.1) 

-'[ 1 / 2 . 1 ) 

then we can express the operator Cq in the simpler form 

Cc,Kf{x) := cxf{x)-\- [ {f{y\x)-f{x))K{dy) (3.2) 

•1 [1/2.1) 

with c = b-\- i)(l~2/)-^(bj/)- This situation was considered in depth by Doumic and 

Escobedo [11], and most of the results of this section should be viewed as extensions of 
those in [11] to the case when either a > 0, or RT fulfills (1.5) but not (3.1). Furthermore, 
the case c < 0 was considered by Haas [13] using the same method that we employ 
below. 

3.1. Main resnlts 

The key observation in Lemma 2.2(i) that power functions hq are eigenfunctions of 
the operator £ 0 , underlies the analysis of the homogeneous case. Specifically, if we 
knew that (/it)t>o solves (1.2) with f = hq for q > 2, then the Mellin transform of ptt, 
Mt{z) = {fj,t,hz-i), would solve the linear equation 


dMt{q -b 1) = K{q)Mt{q -b 1). 


(3.3) 
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Focussing for simplicity on the initial condition /zq = (5i, so that Mo{q) = 1, we would 
find 

Mt{q + l) = exp{tK{q)). (3.4) 

In order to check that (3.4) is indeed the Mellin transform of a positive measure, we 
define, for every w G dom k, a new function by shifting k: 

$(^ (<?):='«(w + g) - k(u;), q>0. 

This is a smooth, convex function, and it has a simple probabilistic interpretation, 
which will play a major role throughout. In this direction, recall first that a Levy 
process is a stochastic process issued from the origin with stationary and independent 
increments and cadlag paths. It is further called spectrally negative if all its jumps are 
negative. If ^ := (C(0)t>o is a spectrally negative Levy process with law P, then for 
alH > 0 and d G R, the Laplace exponent $, given by 

E[exp((7e(l))] = exp(<I>(g)), 

is well defined (and finite) for all g > 0, and satisfies the classical Levy-Khintchin 
formula 


<i>{q) = a.q+l-a'^q'^ + [ (e«“ - 1 + 9 xl{|,j,|<i})T(da;), g > 0, 

^ J ( — 00 , 0 ) 

where a G M, <7 > 0, and T is a measure (the Levy measure) on (—c», 0) such that 


' ( — 00 , 0 ) 


(1 A a;^)T(dx) < oo. 


is in fact the Laplace exponent of a spectrally negative Levy process; see [21, Theorem 

8 . 1 ]. 

Lemma 3.1. Let lo G dom/t. Then: 


(i) the function is the Laplace exponent of a spectrally negative Levy process, 
which we will call fuj = {iu;it))t>o- 

(a) for every t > 0, there exists a unique probability measure on (0, oo) with 
Mellin transform given by 

Mp^f'^q + 1) := [ a;V^'(dx) = exp(t$<^((7)), g > 0. (3.5) 

J (0,oo) 

The family of measures has the representation = P(exp(^i^(t)) G •), fort > 0. 

Proof. We prove both parts simultaneously, and focus first on the case w = 2, where 
we write $ := $ 2 . We can express $ in the form 


<^{q) = aq‘^ + b'q+ [ {y'^ - 1 + q{l - y))y‘^K(dy) + [ ((1 


[ 1 / 2 . 1 ) 


l){l-y)^K{dy) 

(3.6) 
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with 

5 '= 30 + 6 +/ {l-y){l-y‘^)K{dy). 

J [1/2,1) 

Let us denote by A(da:) the image of y'^K{dy) by the map y x = Iny, and n(da;) 
the image of {l — y)‘^K{dy) by the map y ^ x = ln(l —y). Then, thanks to (1.5), A is a 
measure on [— In 2,0) with J a;^A(da:) < oo, and If is a finite measure on (—oo, — In 2], 
and there are the identities 



l + q{l-y))y^K{dy) 


[- ln2,0) 


1 + g(l - e"))A(da;) 


and 


/ {{l-yy-l){l-yfK{dy)=j (e«"-l)n(dx) 

^[1/2,1) ^f-oo,-ln2l 


( —oo, — In 2] 


This shows that $ is given by a Levy-Khintchin formula, and therefore, <I> can be 
viewed as the Laplace exponent of a spectrally negative Levy process ^ = (C(t))t>o 
(see Chapter VI in [1] for background), i.e.. 


E{exp{q^{t))) = exp{td>{q)), t,q>0. 


We conclude that (3.5) does indeed determine a probability measure pt which arises 
as the distribution of exp(^(t)). 

Finally, if w 2, we observe that the function can be written 

‘hi,;)?) = $(y + w — 2) — <I>(a; — 2), which implies that is given by an Esscher 
transform of <&, and hence is also the Laplace exponent of a spectrally negative Levy 
process; see [16, Theorem 3.9] or [21, Theorem 33.1]. 


Remark 3.1. More generally, if C denotes a random time having an exponential 
distribution, say with parameter k > 0, which is further independent of 5, then the 
process 



if t < C 
if t > C 


is referred to as a killed Levy process. Note that if we set = k + d)(( 7 ), then 


]E[exp(y^tW)] = exp(t4>t (<?))> 


with the convention that exp(y^j(t)) = 0 for t > (f. So Lemma 3.1 shows that whenever 
k{uj) < 0, the function q i —>■ k{uj + q) can be viewed as the Laplace exponent of 
a spectrally negative Levy process killed at an independent exponential time with 
parameter —k(w). 

Recall from Lemma 2.1(ii) that 2 £ dom/t always; we will write pt for p). . Since pt 
is guaranteed to exist, this collection of measures will play a particular role in the case 
a = 0. We stress that in the cases a < 0 and a > 0, we will need to choose different 
values of to, and the notation pt will then refer to a different distribution. 

We point out the following property of the probability measures , which essen¬ 
tially rephrases Kolmogorov’s forward equation. 
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Corollary 3.1. The family of probability measures (/9[“^)i>o defined in Lemma 3.1 
depends continuously on the parameter t for the topology of weak convergence. Further, 
for every g G C“(0, oo), the function t i—>■ {p\’\g) is differentiable with derivative 
9t{pt\g) = {p\f\Ag), where 

Ag{x) := Co{hujg){x) - K(uj)g(x), x > 0. 

Proof. Recall that every Levy processes fulfills the Feller property and in partic¬ 
ular, for every function (p G Co(IR.), the map t i—E((/?(5i^(t))) is continuous. Taking 
p>{x) = g{e^) with g G Co(0, oo) yields the weak continuity of the map 1 1 -?- p^f^\ 

Further, it is well-known that the domain of the infinitesimal generator of a Levy 
process contains C“(IR) (see, e.g.. Theorem 31.5 in Sato [21]), and it follows similarly 
that for g G C“(0, oo), the map t i-)- {p^f^\g) is differentiable. To compute the 
derivative, that is to hnd the infinitesimal generator, take q > 0 and recall that 
hq{x) = x’^ for a; > 0. Then simply observe from (3.5) that 

dt{pt,hq) = $^(g)exp(t$^(g)) = (p[‘^',$<^((7)/iq) = {p^f"\K{q + uj)hq-K{uj)hq). 

Using Lemma 2.2(i), we can express 

K{q + Uj)hq = h-cjK{q + Uj)hq+uj = h_,^Co{hujhq), 

which shows that dt{pf^\ hq) = {p\f’\Ahq) for all <? > 0. That the same holds when hq 
is replaced by a function g G now follows from standard arguments, using linear 
combinations of hq. 

We would now like to invoke Lemma 3.1 to invert the Mellin transform (3.4), 
observing that (using (3.5) with p = p'A) 

exp(tK(( 7 )) = exp(t«;(2))(pt,a;‘?“^), 

and conclude that 

pt{dix) = AA‘^''x~'^pt{dx). 

However hq ^ C“(0, oo) and we cannot directly apply this simple argument. Nonethe¬ 
less we claim the following. 

Theorem 3.1. The equation (1.2), for f G C(?°(0,oo) and with C = Cq, has a unique 
solution started from po = Si, given by 

Ptidx) = e*’^^‘^^x~‘^pt{dx), t >0, 

where pt is the probability measure on (0, oo) defined by (3.5) for to = 2. 

Remark 3.2. In particular, the unique solution in Theorem 3.1 fulfills 
{pt,hq) = exp(t«;(g)), as we expected from (3.4). From a probabilistic perspective, 
this does not come as a surprise. In [3], a homogeneous growth-fragmentation stochas¬ 
tic process Z{t) = {Zi{t), Z 2 {t),...) was constructed whose evolution is, informally 
speaking, governed by the stochastic growth-fragmentation dynamics described in 
the introduction. Using a spine technique, it may be shown (we omit the proof) 
that the solution {pt)t>o has the representation {pt,f) = s-iiy 
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/ for which the right-hand side is finite; and in [3, Theorem 1], it is proved that 
= exp(tK(g)) for all q > 2. We offer a more detailed discussion of the 
spine technique in section 6. 

Proof of Theorem 3.1. It is straightforward to check that = e*'^^^^a:“^pt(da;) 

is indeed a solution. Specifically, we deduce from (3.5), that 

(/rt, hq) = exp(tK(2)) exp(t $((7 - 2)) = exp(t«;( 9 )). 

We thus see that {pt)t>o solves (1.2) with C = Cq and f = hq for every <7 > 0, and it 
follows from classical properties of the Mellin transform that this entails that {iit)t>o 
solves (1.2) more generally for all / G C“(0, oo). 

Conversely, given a solution to (1.2) with /io = (5i, set 

ptidx) = e~*'^^'^^x'^fj,tidx). 

Take g G C“(0, oo) and define f{x) = x‘^g{x) for a; > 0, so / G C“(0, oo). Then we 
have {pt,g) = /) and 

dt{pt,g) = -K{2){pt,g) 

that is, 

dt{pt,g) = {pt,Ag), (3.7) 

with 

Ag{x) = x~'^Cof{x) - K{2)g{x), 

as in the notation of Corollary 3.1. We can thus interpret (3.7) as Kolmogorov’s forward 
equation for the infinitesimal generator of the Feller process (exp(^(t)))t>o. This will 
in turn enable us to identify pt = Pt- 

To be precise, examining the proof of [12, Proposition 4.9.18], we see that (3.7) for 
all g G C“(0, oo) has at most one solution (in the sense of a vaguely right-continuous 
collection of measures {pt)t>o) so long as the image of C“(0, oo) by A — ^ is separating 
(see [12, p. 112]) for each A > 0. Since A is the generator of a Feller process and 
C“(0, oo) is a core (cf. Theorem 31.5 in Sato [21]), we know that the image of C“(0, oo) 
by A — Al is a dense subset of Cq, and this implies that it is separating. If {pt)t>o is 
a collection of measures solving (3.7), then for any g G C“, the function t i-G 
is right-continuous. Hence, the solution of (3.7) restricted to is unique, and this 
transfers to (1.2). 

3.2. Some properties of solutions 

We next present some properties of the solution identified in Theorem 3.1, by means 
of translating known results on Levy processes. 

We first point out that, depending on whether (3.1) holds and a = 0, the support 
of the solution pt is bounded or not. Specifically, if a = 0 and (3.1) holds, we set 

d:=b-\- [ (1 - y)K{dy), 

and otherwise d = oo. It is easy to verify that d = limg_>oo q~^K{q). 
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Corollary 3.2. If a = 0 and (3.1) holds, then for every t > 0, is the supremum of 
the support of p,t, i.e., we have for every e > 0, 

cx))) = 0 and - e, e*'^]) > 0. 

Proof. The spectrally negative Levy process ^ = ^ 2 , arising in Lemma 3.1, has 
bounded variation with drift coefhcient d exactly when the conditions of the result hold, 
and it is then well-known that td is the maximum of the support of the distribution 
of £,{t). Therefore we have pt((e*'^, 00 )) = 0 and — e, > 0, and our claim 

follows from Theorem 3.1. 


In the case when the assumptions of Corollary 3.2 are not fulfilled, we have the 
following large deviations estimates for the tail fit{x) ■= of pt. Recall that 

K is a convex function, and observe that limg_>+oo K!{q) = -l-oo when either a > 0 or 
(3.1) fails. Thus for every r sufhciently large, the equation k!{ q) = r has a unique 
solution which we denote by 0{r), and the Legendre-Fenchel transform of n is given by 


K*{r) := supjrg — K{q)} = r9{r) — K,{6{r)). 

q>0 


Corollary 3.3. Suppose that a > 0 or (3.1) fails. Then for every r > 0 sufficiently 
large, we have 

lim t~^\Yifit{e^'^) = —K*{r). 

t—¥00 

Proof. This follows easily from the identity (ptffiq) = exp(t«;(( 7 )) by adapting the 
classical arguments of Cramer and Chernoff; see, for instance. Theorem 1 in Biggins 

[7]. 


The estimate of Corollary 3.3 can easily be reinforced by using the local central limit 
theorem. Here is a typical example (compare with Theorem 1.3 in [11]). 

Corollary 3.4. Suppose that a> 0 or (3.1) fails, and further that K'{q) < 0 for some 
q. Then 6(0) is well-defined, 0 < k"(—0(0)) < 00 , and for every f G Cc, we have 


(TtJ) ~ 


giK(e(o)) 

^/2'KtK"{0{O)) 





Proof. The first assertion about the existence of 0(0) and «;"(0(O)) are immediate 
from the convexity of k and the fact that lim +00 k = -l-oo. 

The function l>(g) := K{q 0(0)) — k( 0(O)) = $(g -I- 0(0)) — $(0(0)) is the Laplace 
exponent of a spectrally negative Levy process (5(t))t>o which is centered and has 
finite variance «;"(0(O)). Further, we see from Esscher transform and Theorem 3.1 that 

^,(dx) = e‘«(^(°»x-^(°)p(exp(e(t)) e dx). 


Our claim then follows readily from the local central limit theorem for the Levy process. 

Corollary 3.5. If a > 0 or the absolutely continuous component of K{dy) has an 
infinite total mass, then pt is absolutely continuous for every t > 0. 

Proof. Using Sato [21, Theorem 27.7 and Lemma 27.1], it follows from the assump¬ 
tions of the statement that the one-dimensional distributions of the Levy process f{t) 
are absolutely continuous for every t > 0. Our claim follows from the representation 
in Theorem 3.1. 
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4. The self-similar case 

We now turn our attention to the growth-fragmentation equation (1.2) for C — Ca 
given by (1.4) and a ^ 0. We first point out that the function k is non-increasing if 
and only if a = 0, the dislocation measure K fulfills (3.1), and 

b+ f (1 -y)K{dy) < 0. 

“'[ 1 / 2 . 1 ) 

In this case, the operator can be expressed in the form (1.3) with c < 0, and (1.2) 
is then a pure fragmentation equation as studied by Haas [13]. To avoid duplication 
of existing literature, this case will be implicitly excluded hereafter. 

Recall the notation hq{x) = x® for a: > 0. In the self-similar case, power functions 
are no longer eigenfunctions of the operator however, there is the simple relation 

Cahq = K{q)hq+a, G dom k; (4.1) 

see Lemma 2.2(i). Hence, if (1.2) applies to power functions, the linear equation (3.3) 
for the Mellin transform Mt{z) = {nt, hz-i) in the homogeneous case has to be replaced 
by the system 

q) = K{q)Mt{l-^ q-^ a). (4.2) 

We make the fundamental assumption, that 

inf K{q) < 0, (4.3) 

q>0 

which is implicitly enforced throughout this section. The role and the importance 
of (4.3) shall become clear in the sequel. Recall that k is a convex function on K, 
and is ultimately increasing, since we are excluding the case when k is non-increasing 
throughout section 4. Hence, condition (4.3) ensures the existence of a unique G R 
with 

k{oj^) = 0 and > 0. (4.4) 

We refer to as the Malthusian parameter. 

The sign of the scaling parameter a plays a crucial role, and we shall study the two 
cases separately, even though some ideas are similar. 

4.1. The case a < 0 

We now focus on the case a < 0. We start by observing that the existence of a 
Malthusian parameter enables us to view (4.2) as a closed system for an arithmetic 
sequence, and thus to solve it. 

Lemma 4.1. Consider a sequence of functions M,(l -j- q): [0, oo) (0, c») for 
q = Lo^ — ka, k = —1,0,1, ..., with Mq{ 1 q) = 1. Suppose that (f.2) and (4-3) hold 
and recall that uj+ is the Malthusian parameter defined by (4-4)- Then Mt{l + w^.) = 1 
for all t > 0 and for fc = 1, 2 ,..., we have 

M,(1+-1„) = 1 + y: "K-»t)-»K-«(t-g+i)) 

e=i 
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Proof. The equation (4.2) applied to the Malthusian exponent implies that the 
function t i-T Mt{l + is constant. We can then solve (4.2) for q = — ak and 

fc = 1, 2,... by induction in order to obtain the given formula. 

In comparison with the homogeneous case, Lemma 4.1 is a much weaker result than 
(3.4), as we are not able to compute the whole Mellin transform of a solution, but merely 
its moments for orders forming an arithmetic sequence. There is hence an additional 
crucial issue: it does not suffice to find a family of measures having the desired 
moments, but also to ensure that the moment problem is determining. It turns out 
that moment calculations which were performed in [6] for self-similar Markov processes 
enable us to solve the moment problem in Lemma 4.1, and check that this indeed yields 
a solution to (1-4). Similar calculations also point at a rather surprising result, namely 
that the self-similar growth-fragmentation permits spontaneous generation! 

Theorem 4.1. Assume (4-3) and a < 0. 

(i) For every t > 0, there exists a unique measure on (0, oo) such that 
{/If, ha,^) = 1 and for every integer k >1, 


{Ft,hu 


-ka 


^ ^ ^ k{uj^ — ak) ■ • ■ — a(k — £ + 1)) 


£l 


In particular, /Tq = di and the family it^t)t>o solves (1.2) for all f € C“(0, oo) 
when C = Col is given by (1.4)- 

(a) For every < > 0, there exists a unique measure on (0, oo) such that (y^, h ^^^) = 1 
and 


(To f^0J,—ka) — t 


i{uj^ —«)••• k(uj^ — ak) 
k\ 


for every integer k > 1. 


If we further set 7 q = 0, then the family (jf)t>o solves (1.2) for all f G C“(0, oo) 
when C — Col is given by (1-4)■ 

Proof. (i) Let us define 4)^ := = K(--|-a;+), which, as we saw in Lemma 3.1, is 

the Laplace exponent of the Levy process := ■ Observe that 

<1>'^(0) = k'(w^) > 0, so this Levy process has a strictly positive and finite first 
moment. 

Proposition 1 in [6] then ensures, for every t > 0, the existence of a unique 
probability measure on (0, oo), such that for every integer A: > 1, 


(Pt> ^ — ak) 


1-f 


k{uj. 




ak) ■ ■ ■ k{uj^ 

w 


a{k-£+!)) f 

-A , 


so that in particular p^ = d\. Thus we may set Pt{dx) = x~‘^+ p^)dx), and then 
{Pf, hcj^-ko:) is given as in the statement for every integer k > 0. That this 
determines derives from the uniqueness of p^. 

Now, using (4.1), we immediately check that {pt, satisfies (4.2). It then 

follows that {pt)t>o solves (1.2) for every / G (0, oo) (recall that the probability 
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measure is determined by its entire moments). Finally, the map 1 1—1 {p^, h-ak) 
is continuous, and we deduce that {pt)t>o is vaguely continuous (using again the 
fact that Pf. is determined by its moments (p^, h-ak) for k G N). Hence the same 
holds for (pt)t>o- 

(ii) Recall from above that = k{uj^ + •) is the Laplace exponent of a spectrally 
negative Levy process which has strictly positive and finite first moments. Propo¬ 
sition 1 in [ 6 ] ensures for every t > 0 the existence of a unique probability measure 
TTj on ( 0 , oo) such that its moments are given by 

(7r„ h-ak) = t - - - = t - - -, 

for k = 1,2,..., and this determines tt^. It follows immediately that ( 7 rj)t>o is 
vaguely continuous (recall that ttq = 0 ). 

We then define for t > 0 


7 j(dx) = a; ‘^+ 7 r((dx), a; > 0 , 


so 


(dt) huj —ka) — t 


i{uj^ —«)••■ k(uj^ — ak) 

k\ ■ 


Then (4.1) entails that for every integer fc > 1, there is the identity 

,k-i'^i^+ -a)---K(w+ - ak) 


(bi5 hoj —ka) — i 


(fc-1)! 

= k{ui^ - afc)(7t>^cj+-(fc-l)a) = {lt,t^h^^-ka), 


and the conclusion follows just as in (i). 

Theorem 4.1(ii) entails that uniqueness of the solution fails when one only requires 
(1.4) to be fulfilled for all / G C“(0, oo), which contrasts sharply with the results of 
Haas [13] for the pure-fragmentation equation. We conjecture that the solution p,^ 
given in Theorem 4.1(i) is minimal, in the sense that if {p,t)t>o is another solution with 
the same initial condition = 5i, then p^ < pt for every t > 0. We also stress that 
uniqueness of the solution can be restored by requiring ( 1 . 2 ) to hold for the functions 
hq with q > a] see Lemma 4.1 and Theorem 4.1(i). 

We now present a different approach to Theorem 4.1. In the homogeneous case a = 0, 
we saw in the preceding section that the equation (1.4) bears a close relationship with 
certain exponential Levy processes. It turns out that in the self-similar case with a < 0, 
the vital connection is with positive self-similar Markov processes, and is made via the 
Lamperti transformation which associates these with the class of Levy processes. We 
first provide some background in this area. 

A positive self-similar Markov process {pssMp) with self-similarity index 7 G M 
is a standard Markov process R = {Rt)t>o with associated filtration {^t)t>o and 
probability laws (Pa;)a;G( 0 , 00)1 on [ 0 , 00 ], which has 0 and 00 as absorbing states and 
which satisfies the scaling property, that for every x,c> 0, 


the law of {cRtc-<^)t>o under P^, is Pea;- 
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Here, we mean “standard” in the sense of [ 8 ], which is to say, is a complete, 

right-continuous filtration, and R has cadlag paths and is strong Markov and quasi- 
left-continuous. 

In the seminal paper [17], Lamperti describes a one-to-one correspondence between 
pssMps and (possibly killed) Levy processes, which we now outline. It may be worth 
noting that we have presented a slightly different definition of pssMp from Lamperti; 
for the connection, see [ 22 , § 0 ]. 

Let S{t) = fg (Ru)~'’'du. This process is continuous and strictly increasing until R 
reaches zero. Let (T(s))s>o be its inverse, and define 

ys= log Rt(s) s > 0 . 

Then rj := (jis)s>o is a (possibly killed) Levy process started at position logx, possibly 
killed at an independent exponential time; the law of the Levy process and the rate 
of killing do not depend on the value of x. The real-valued process rj with probability 
laws (Py)ygR is called the Levy process associated to R, or the Lamperti transform of 
R. 

An equivalent definition of S and T, in terms of rj instead of i?, is given by taking 
T(s) = fg exp(jr/u) du and S as its inverse. Then, 

Rt = exp(r/s(t)) 

for all t > 0, and this shows that the Lamperti transform is a bijection. A useful 
fact is that, as a consequence of the definitions we have just given, it holds that 
dt = exp(-jr/g(t))dS(t). 

Most of the literature on pssMps (including Lamperti’s original paper) assumes that 
7 > 0, and much of it is also given for 7 = 1 . Indeed it is easy to change the index 
of self-similarity. If i? is a pssMp of index 7 and corresponding to the Levy process 
r], then, for any 7 ' G M, the process i?’’' = {R] )t>o is a pssMp with index 7 / 7 ', 
corresponding to the Levy process y'ry. It is also useful to note that the time-changes 
appearing in the Lamperti transformation are a.s. equal for R and i?'*' . We should 
point out that the case 7 = 0 is special, since in this case the time-change does not 
have any effect, and the pssMps of index 0 are just exponential Levy processes. 

Note that, if the Levy process process p is killed at time f, then we define Rt = 0 
for t > T{C,) if 7 > 0, and Rt = -too for t > T{f) if 7 < 0. 

Recall Lemma 3.1, and define := <I>(^ = k(- -I-w^), which is the Laplace exponent 
of the spectrally negative Levy process := Let us denote by the pssMp 
with index —a associated to by the Lamperti transformation. Note that, because 
^_i_ has positive mean, the process never reaches the absorbing boundaries 0 or 
-boo. We define the measure to be the distribution of X^{t) under Pi, that is, the 
probability measure on ( 0 , 00 ) defined by 

(p„/) = Ei(/(Xjt))), /GCo(0,(^) 

and give first the following analogue of Corollary 3.1: 

Lemma 4.2. The family of probability measures {pt)t>o depends continuously on the 
parameter t for the topology of weak convergence. Further, for every g € C“(0, 00 ), the 
function 1 1 —>■ (Pt,g) is differentiable with derivative dt{pt,g) = where 

A^^^g{x) := x~‘^+Ca{hu:^g){x) , 


a; > 0. 
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Proof. The first assertion follows easily from the Feller property of self-similar 
Markov processes; see Theorem 2.1 in Lamperti [17] and the remark on page 212. 
In order to establish the second, we work with the Levy process The exponential 
Levy process exp(^^(-)) is a Feller process in (0,oo), and the same calculation as in 
the proof of Corollary 3.1 shows that its infinitesimal generator is given by 

■^+9{x) = x~‘^+Co{haj^g)ix), g G C“(0,oo). 

According to Dynkin’s formula (see, e.g., Proposition 4.1.7 of [12]), for every 
g € C“(0, oo), the process 

5(exp(C+(t))) - f A_,5(exp(C+(s)))ds 

Jo 

is a martingale. Recall that by definition, arises as the transform of exp(^^(-)) 
by the time substitution S', which is given as the inverse of the additive functional 
fo (®xp(^+(s)) ds, and we have the identity 

rS(t) ft 

ff(^+W)-/ -4+ff(exp(C+(s)))ds = 5 (X+(t)) - / /i„(X_^(s))A_,g(X_^(s))ds. 

Jo Jo 

A priori, the time-substitution above changes a martingale into a local martingale. 
However, using Lemma 2.1(iii) and the fact that a < 0, we see that g := haA^g 
is bounded, and it follows that the process 

g{X^{t))- f A[^'>g{X4s))ds 
Jo 

is a true martingale. Taking expectations, we arrive at 

{Pt,9)-9i0= [ (Ps,T“^5)ds, 

Jo 

and our claim follows. 

The connection with solutions of the growth-fragmentation equation is the following: 


Corollary 4.1. Let 

jlt = h-u:^Pt, t>0- 

Then, (fif)t>o is equal to the solution {g,^)t>o of the growth-fragmentation equation 
appearing in Theorem 4-Hi). 

Proof. We deduce immediately from Lemma 4.2 that for every / G 

dtiftJ) =9t{Pt,hZlf) = {PtHfl^af) = (PtHaf), 

that is, the family {ft)t>o (1-2) with C — La- That coincides with the measure 

appearing in Theorem 4.1(i), and that the notation for the distribution of X^{t) is 
consistent with that used in the proof of Theorem 4.1(i), follows from Proposition 1 of 
[ 6 ]. 
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This approach could also be adapted to prove the existence of (jf)t>o using the 
process started from zero, and indeed, this will be our method for the case 

a > 0 in section 4.2. 

We conclude the section by offering some results on the asymptotic behaviour of the 
solution (/it)t>o given by the previous theorem. 

Our first result in this direction indicates that, thanks to the self-similarity property 
(2.1) of the equation (1.2), the solution starting from zero given above can be used to 
describe the asymptotic behaviour of as t —?■ oo. 

Proposition 4.1. For any / G Cb(0, cxj), 

J /(t-i/l“lz)z-+Mt(d2) ^ J /(z)z“+7i(d^). 

Proof. Since k'{oj+) > 0, it is possible to extend the definition of in order to 
allow it to start from (0) = 0, such that it is a Feller process on the state space 
[0, oo); this is a consequence of [5, Theorem 1]. For a; > 0, we will denote by Px the 
law of the process with X^ (0) = x. 

It then follows from the scaling property that Ex[f{t^^°‘X_^{t))] = E,j,(i/c [/(A"^ (1))], 
and then the convergence 

E,[/(ti/“X^(t))] ^ Eo[/(Xjl))], t ^ oo, 

follows from the scaling property of X^. 

Finally, we know from the reference [6], which we used in the proof of Theorem 4.1, 
that the measures x‘^+ii^{dx) and a;‘^+ 7 ((da:) are, respectively, equal to Pi(X+(t) G dx) 
and Po(^+(t) G dx). Our claim follows immediately. 

We remark that the statement of the proposition can easily be extended to solutions 
based on (/r^) whose initial value is a measure with compact support in (0, oo). 

Suppose now that the equation K{q) = 0 has two solutions, uj_ and u;+, with uj_ < 
Then we can say a little more. Let X_ be the (—a)-pssMp associated with the Levy 
process := having Laplace exponent <I)_ := . Recall that we say 

the Levy process is lattice if, for some r G M, the support of ^-(1) a.s. lies in rZ; 
otherwise, we say that is non-lattice. If we define the random variable 

poo 

1 = 

^0 

then it is known from [19, Lemma 4] that, so long as is non-lattice, 

lim t(‘^+-‘^-)/l“IPo(/ >t)=C, 

t—>oo 

for some 0 < C < oo. We obtain from Haas and Rivero [14] the following result. 
Proposition 4.2. Let / G Co(0, oo) and assume that is non-lattice. Then, 


J x‘^-f{t ixl“l)/xt(dx) 
MI > t) 


f{x)iy{dx), 


t oo. 


where v is the distribution of the random variable )/\a\ equation (13) of [1)]. 
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Proof. As remarked in the proof of Proposition 4.1, the equation 

x'^+fj'tid^) = Pi(A+(f) G dx) 

holds as an identity of probability measures, where A_|_ is the (—a)-pssMp correspond¬ 
ing to the Levy process with Laplace exponent $+ = . We now wish to use the 

‘Esscher transform’ for pssMps, which is essentially obtained by standard arguments 
from the Esscher transform of Levy processes given in [16, Theorem 3.9] (see, for 
instance, the discussion around [9, Theorem 14] for an application in the context of 
pssMps.) This allows us to perform a change of measure to switch from the process 
A+, related to the Laplace exponent $+, to the process X_, related to the Laplace 
exponent 4>_ = <!)+(• -\-oj_ — uj+). Specifically, we have: 

=x^--^+x‘^+Ht{dx) =x^--‘^+PiiX+{t) G dx) = Pi(A_(t) G dx), 

for X > 0. The process (which corresponds to A_) is a Levy process with non¬ 
monotone paths and which satisfies the conditions of Haas and Rivero [14, Theo¬ 
rem 1.6]. Applying this theorem gives the result. 

4.2. The case a > 0 

In the case a > 0, the equation (4.2) for the Mellin transform is unfortunately much 
less useful, for the following reasons. Firstly, the analogue of Lemma 4.1 would require 
to assume that {/it^hq) < oo for all q sufficiently negative. Roughly speaking, this 
would force the scarcity of small particles, and this phenomenon only occurs for a very 
restricted class of dislocation measures K (informally, dislocations should not generate 
too many small particles, and in particular the total intensity of dislocations must be 
finite). Secondly, even if one were able to get an expression for the (negative) moments 
(/ij, huj^-ka) with fc G N, this moment problem would be in general indeterminate, and 
the arguments used in the preceding section would thus collapse. 

Nonetheless, we have just seen from Lemma 4.2 that, for a < 0, self-similar growth- 
fragmentation equations have a close connection with certain self-similar Markov pro¬ 
cesses, and using the intuition that we gained, we are able to offer a very similar set of 
results when a > 0. Recall that k: [0, oo) (—00,00] is a convex function, and that, 
since we are assuming (4.3), we may pick w > 0 such that k := —k(w) > 0. Then, the 
function 

$1(9) := -kk = K(w-k(7), g > 0, 

is the Laplace exponent of a spectrally negative Levy process killed at an independent 
exponential time with parameter k, say and we denote by A| the (—a)-pssMp 
associated with via the Lamperti transformation. Note that A| hits the absorbing 
state - 1-00 by a jump. 

We write for the sub-probability measure on (0, 00 ) induced by the distribution 
of A|(t) and study its properties in the following result, which mirrors Lemma 4.2. 

Lemma 4.3. Suppose (4-3) holds and a > 0. Then we have, in the notation above: 

(i) E(sup(>Q X|(t)‘?) < 00 for all 0 < q < — uj. 
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(in) The family {p^)t>o depends continuously on the parameter t for the topology 
of weak convergence. For every g € C“(0, oo), the function t i—>■ {pt,g) is 
differentiable with derivative dt{pt,g) = {pf,A^^°‘'^g), where 

'■= x~‘^Ca{hu,g){x), a: > 0. 


(iv) {pf)t>o solves the above equation also for g{x) = x® with 0 < q < uj^ — uj. 

Proof. (i) From the very construction of the overall supremum 

Xf := supt>oXt(^) is given by exp(^t). with := supt>oCt(^)- We 

infer from Corollary VIL2 in [1] that has the exponential distribution with 
parameter — w (which is the positive root to the equation ^.\{q) = 0), and our 
claim follows. 


(ii) We begin with the following calculation, using the discussion of pssMps in the 
preceding section. Recall that S is the time-change appearing in the Lamperti 
transform relating X.f and and that there is the identity dS{t) = exp(Q;^l(t))dt. 
We have therefore 


POO POO 

/ X^{ufdu= / eP«t(^(“»du 

Jo Jo 

pOO pOO 

= e(P"“)«t(S(“)) dS'(M) = / 

Jo Jo 


But now we can consider the expectation: 


E 


X 



X^{uy du 


xP-°E 


j(p-a)^As) ds 


xP-a 

k{p — a + w)' 
oo, 


if k(p — a + ui) < 0, 
otherwise. 


We complete the proof by recalling the definition of . 

(hi) Just as in the proof of Lemma 4.2, the first assertion follows from the Feller 
property of self-similar Markov processes, and the process 


Nt := g{Xyt)) - 5(1) - f M;“^5(^t(s))ds (4.5) 

is a local martingale for every g G C^(0, oo). We will show that 
lE(sup^<t|fVs|) < oo, t>0, 

which implies that TV is a true martingale; see [18, Theorem 1.51]. 

Since g is bounded, certainly sup 5 <(| 5 (Ar.|.(s)) — 5 (a;)| is in L^(P). In constrast 
to Lemma 4.2, the function g may be unbounded for a > 0; however, we do 
know from Lemma 2.1(iii) that is zero on some neighborhood of 0 and, for 

any q G doms:, fulfills g = o(a;^''‘“““) as a; —)■ oo. 
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We let w < <7 < and keep it fixed for the rest of the proof. For some K > 0 
we then have 


E 


sup 

U<t 




< E 




(4.6) 


< t sup + E 

[0,K] 


Jo 


Setting p = q a — u) in part (ii), we see that the right-hand side is finite. 

We have thus shown that is a true martingale, and 

E (/ < oo- 

Taking expectations in (4.5) and applying Fubini’s theorem, we obtain 

(dt,5)-ff(l) = / 

Jo 

which completes the proof. 

(iv) This part is proved by setting g(x) = x® in the previous part, as follows. Using 
the Markov property one sees immediately that the process 

Mt = e«ft(‘) - 1 - k{uj_ +q) f ds, t > 0, 

Jo 

is a martingale in the filtration of for every g > 0. Applying the same reasoning 
with the time-change as in Lemma 4.2, it follows that 

Nt = X^{ty - \ - k{uj + q) [ A:t(s)«+“ds, t>0, 

Jo 

is a local martingale. (For our choice of g, we have A^'^^g^x) = k{uj q)x'^^°‘, so 
this is consistent with the proof of part (hi).) We observe that 

poo 

sup|Ait| < 1-I-sup A'.|.(t)^ — k(w- f g) / X|(s)'^+“ds. 

t>o t>o Jo 


We now apply directly parts (i) and (ii) of this lemma in order to show that 
E[sup(>Q|Alt|] < oo. This is a sufficient criterion for to be a uniformly integrable 
martingale (see [18, Theorem 1.51]), and the remainder of the proof follows in the 
same way as in part (hi). 

We can then repeat the calculations which were made after the proof of Lemma 4.2, 
and arrive at: 

Corollary 4.2. Suppose that (4-3) holds. Define = h-c^p^., for t > 0. Then 
the vaguely continuous family of measures {p,f)t>o solves (1.2) with C = Ca both for 
f G C“(0, c») and for f = hq, for any to < q < LO+. 
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An interesting contrast with the case a < 0 is that we do not show that solves 
(1.2) for all power functions. 

We now give the basis of a solution to the growth-fragmentation equation starting from 
the zero measure; in this case, the solution should be interpreted, not as spontaneous 
generation from infinitely small masses, but as starting from infinite mass and breaking 
apart instantaneously. Recall that the equation K{q) = 0 has at most two solutions. 
More precisely, we have already seen that there is always a unique solution with 
> 0 (this is the Malthusian exponent defined by (4.4)). When a second solution, 
say w_, exists, then w_ < and k'{ui_) G [— oo, 0 ). We give the results under the 
assumptions: 

the equation K{q) = 0 with q >0 has two solutions uj_ < and k'{uj_) > —oo 

(4.7) 

which is thus stronger than (4.4). We write for the spectrally negative Levy process 
with Laplace exponent 4>_(g) := = K{q + uj_), and then X- for the pssMp with 

index —a associated to by Lamperti’s transform. 

Lemma 4.4. Assume that (4-V holds. Then there exists a cddldg process {X{t))t>o 
with values in (0, oo) and limt^o-i- A!{t) = oo a.s. such that: 

• For every s > 0, conditionally on X{s) = x, the shifted process {X{s + t))t>o 
has the law Px of the pssMp X_ started from x. 

• For all 0 < e < (a;+ — to_)/a, there is c{e) < oo such that 

=c{e)F-\ t>0. 

Proof. Let Y denote the pssMp with self-similarity index a associated to the Levy 
process —f-, so Y has the same law as 1/X-. Because 

E(-^_(l)) = -4)'_(0-|-) := m G (0, oo), 

[5] shows that 0-1- is an entrance boundary for Y, i.e., there exists a cadlag process 
(3^(0)t>o with values in (0, oo) and limt_,,o+3^(t) = 0 a.s., such that, for every s > 0, 
conditionally on 3^(s) = y, the shifted process iyis + t))t>o has the law Y started from 
y. Our first claim follows by setting X{t) = l/y{t). 

Further, according to Theorem 1 in [5], there is the identity 

V / am 

where I := exp(a^_ (5))d5. It thus follows 

=c{e)F-^ 

where c{e) = E (P) /am G (0, oo]. 

For every 0 < e < (u;+ — uj__)/a, the Laplace exponent q ^ ^-{aq) of af- fulfills 
$_(«£) < 0, and according to Lemma 3 in Rivero [20], this ensures that E (/®) < oo. 

We next further require that u!_ G (dom/t)°. This is only a little stronger than the 
condition k'(w_ ) > —oo, which is necessary and sufficient for X to exist. We write 
for the operator given in Lemma 4.3(iii) for tv = uj_, and deduce the following. 
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Corollary 4.3. Assume that (4-7) holds and that G (domK)°. For t > 0, write ivt 
for the distribution of X{t). Then, for every f G C“(0,oo), we have 

ds < oo, 

and there is the identity 

{'^t:f)= [ {TTs,A^°^f)ds. 

Jo 

Proof. Recall from Lemma 4.3(iii) that huj_A^°'^f = Ca(huj_f), so Lemma 2.1(iii) 
and the assumption that u;_ — «£ G dom k for some e > 0 entail 

It now follows from the preceding lemma that \(tTs, f)\ < C{s‘^~^ + 1), where C is 
a finite constant depend ending only on / and s. Our first claim follows. 

Recall then from the proof of Lemma 4.3(iii) that 

fiX_{t))- f A^^^f{X_{s))ds, t>0 

Jo 

is a local martingale, and thus, thanks to Lemma 4.4, so is 

f{X{t + s)) - fiX{s)) - [ A^°^f{X{r + s))dr , t > 0 
for all s > 0. Since 

E^y \A^°^ f{X{rs))\dr'^ = 1(71^, ^L“V)|dr < oo, 

the above process is actually a true martingale, and taking expectations, we arrive at 
the identity 

(T^tjf) - = y {■Kr,A‘ff'’f)dr. 

With the observation above, we can let s 0+, and since (tts,/) 0 thanks to 

Lemma 4.4, we conclude that 

= f {TTr,A^ff''f)dr. 

Jo 

We conclude this section with the following corollary, which demonstrates the ex¬ 
istence of a solution to the growth-fragmentation equation started from zero mass. 

Corollary 4.4. Suppose that the hypotheses of Corollary 4-3 are fulfilled. Let 
7 t(dx) = x~‘^-7rt(dx) for t > 0, and set 7 q = 0. Then, the family (7j)t>o solves 
(1.2) for all f G C“(0, oo) when C = Ca is given by (1.4). 
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5. Explosion of the stochastic model 

In this section, we discuss the behaviour of the stochastic growth-fragmentation 
process in a simplified setting. In particular, we point out that, when the Malthusian 
hypotheses from section 4 do not hold, this stochastic model may experience explosion, 
in the sense that some arbitrarily small compact sets contain infinitely many particles 
after a finite time. We will focus on the case where a < 0, though similar arguments 
can be made for a > 0. 

Assume that the measure K is a probability measure on [1/2,1) which is not equal 
to di, and denote by F a random variable with law K. Choose c € R such that 
E[log(l - F)] -t c < 0 < E[log F] -t c. 

We now set up the stochastic model. Let lA = lJn>o{-^> where {L, i?}° = {0}. 
We view this as a binary tree, as follows. The root node 0 gives rise to child nodes 
L and R; the former then has children LL and LR, while the latter has children RL 
and RR, and so on. We introduce the ancestry relationship by saying that, for 
individuals u,u' € U, u ^ u' if and only if there exists u" € U such that uu" = u'\ we 
also define the strict relation u ^ u' to mean that u ^ u' but u ^ u'. 

Let V = {L, R}^. This is the set of infinite lines of descent, or rays, in U. For 
instance, LLRRLRLRRRRL ■ ■ ■ G V traces a line of descent starting at individual 0, 
and proceeding to L, then LL, then LLR, and so on. ItuGU and u G V, we say (by 
slight abuse of notation) that u A u if and only if there exists v' € V such that uv' = v. 

To each u G U, we assign, independently of everything else, a lifetime T„ which 
is has an exponential distribution of rate 1, and an offspring distribution F„ which is 
distributed with law K. We then recursively assign positions fu to the individuals in 
U. The root is positioned at a given point x G R, that is, ^0 = x. Its descendents are 
positioned as follows: 

CuL = C« + log(l - Fu)-b cTti and CuR = Cu + log{Yu) + cTu, uGU. 

This gives a model in which each individual lives an exponential time, dies, and (on 
average) scatters one child to the left and one to the right. 

It will be convenient to introduce a model in which individuals also move contin¬ 
uously, as follows. For u G U, define its birth time a„ = death 

time bu = + Tw Its position between those times is then given by 

^uit) = Cti + c(t — Uu), for Uu <t < bu- By another abuse of notation, let us define also 
the positions of a ray: for v GV, let ^vit) = ^uit), where u GU is the unique individual 
with Ou Y t < bu and u -< v. We may now see the model as containing individuals 
which move to the right at constant rate c, until an exponential clock rings and the 
individual dies, scattering offspring to the left. 

The model may also be viewed as a stochastic process 3^ = (3’(t))t>o, with 

yi^) ~ 'y '^exp(^„(t)) l{a„<t<6„}) 

uGU 

taking values in the space JV of locally finite point measures and such that 3^(0) = 6e^. 
With this perspective, it is an important and useful fact that the process has the 
branching property. Loosely speaking, this means that the behaviour of {y{t + s))s>o 
is given by collecting the atoms of y{t) and running from each one an independent 
copy of 3^; for a precise statement and proof see, for instance, [3, Proposition 2]. 
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The process y we have just described is a stochastic model corresponding to the ho¬ 
mogeneous fragmentation equation. In particular, if we define a collection of measures 
ipt)t>o via 


{f,uf)=myit),f)]=E 


/(exp(C«W))l{a„<t<b„ 


'-uGU 


/GC= 


then we obtain a solution to (1.2) with a = 0 and C given as in (1.3); we give some 
more details on this in section 6. This corresponds to the function k satisfying, for 
<7 > 0 , 

K{q)=cq-\-f [y'^-\-{l-yy-i\K{dy) > q(c-\- [ logy K{dy)) > 0. 

Here, the first inequality holds because it holds for the integrands, and the second 
inequality is by our assumption about c at the beginning of the section. Also, k( 0) = 1. 
Thus, infq>o K(q) > 0, and so the Malthusian hypothesis (4.4), which was an important 
assumption in section 4, is not satished for our model. 

We now give the model corresponding to the self-similar equation. To do this, we 
should hrst introduce the notion of a stopping line. We say that S = {Sv)vev is a 
stopping line if: 

(i) For every u G V, S'.!, is a stopping time for the natural filtration of 

(ii) If u GU and v,v' GV such that u v, v', then P(5'^ = Sv> | a„ < 5”^ < 6„) = 1. 
Now, for V € V and a € R, let 


n{s)= [ dr, 

Jo 

and denote its inverse by Sy. Then, S{t) = {Sy{t))y^v is a stopping line for every 
t > 0. li u G hi IS an individual such that, for some v GV with u ^ v, Uy < Sy{t) < 
holds, then we define Suit) to be equal to Syit); by property (ii) of the definition of a 
stopping line, this does not depend on the choice of v. We define 

Xy{t) = exp(^^(S'^(t))), V gV, t > 0, 

and 

J^iji ^ ^ ^exp^,i(5u(t)) ^{au<>S'u(t)<bu}’ t ^ Oj (^■^) 

u^U 

where the sum is over only those u for which S'„(t) is defined. The process X is called 
the a-self-similar fragmentation process. The stopping-line nature of S means that the 
process X retains the branching property; however, it is not clear that it should be 
locally hnite, and indeed, our main result in this section is that X a.s. does not remain 
locally hnite for all time: 

Proposition 5.1. Let a < 0. For any a > 0, there exists some random time a such 
that there are infinitely many individuals of X(a) in the compact set [1,1 -I- a]. 



24 


JEAN BERTOIN AND ALEXANDER R. WATSON 


Proof. We will study rays pk = L^R°° which follow the left-hand offspring for k 
steps, and the right-hand offspring thereafter. Our first remark is that, if we dehne 

tq{v) = inf{t > 0 : Xy{t) =0}, u € V, 


then we have 

pOO 

To := ro(L'^) = Tl<^{oo) = / 

^0 

Since is a Levy process with negative mean (by our assumption on c) we obtain 
that To < oo almost surely. 

Consequently, for any p > 0 there exists some infinite set C such that, for k € C, 
Tpk{^L*‘) G ('^0 “ and Xpi_{Tp^{b]^k)) —>■ 0 as fc —>■ oo. We define 

(v) = sup{t > 0 : Ar„(t) < 1}, u G V, 

which is the last passage time of the level 1 by the process Xy; then, for k G C, we 
have 

LtiPk)=Tp,ibp^k)+LtiRn, 

where L^{R°°) is the last passage time of the level 1 by Xro^, computed for an 
independent self-similar fragmentation process started at Xpi^{b]^k). 

We are therefore reduced to studying hrst passage times of the (—a)-pssMp Xj^oc 
corresponding to a spectrally negative Levy process started from a level x < 0 
with Laplace exponent 


^R°-iq)=cq+[ {y'^-l)K{dy). 

Jih^) 

We seek t such that, with positive probability (not depending on x), L^{R°°) < t. 

The first observation is that (0-I-) > 0, which implies that the process Xf^oc 
can be extended to start at zero; it is then Feller on [0, c»). Furthermore, the pssMp 
drifts to - 1-00 as t —>■ oo. Hence, we pick e > 0 and t > 0 such that Po(L('' < t) > 2e. 
By the Portmanteau theorem, liminfa;_>o Pa;(L('’ < t) > 2s also; and therefore for x 
sufficiently close to zero, Px{Li < t) > e. Applying the Borel-Cantelli lemma to the 
paths referred to above, there are then inhnitely many k G C such that 

Lt (Pk) < Tp^i^Lk) + t <To + t 


with probability 1. 

We therefore have infinitely many paths whose last passage time of 1 occurs in a 
(random) compact interval. In particular, there must exist some finite random time a 
such that for every (5 > 0, there are infinitely many paths which cross 1 for the last 
time in (cr — 5, a). Since the particles are large, the time-change is bounded, in that 
Spk{Lt{Pk) + u) — Sp^{L^{pk)) < u for all m > 0. Furthermore, the processes can 
grow at most at rate c; thus, picking 5 < c~^ log(l -I- a) ensures that, at time a, all the 
selected particles are in [1,1 -f a], which completes the proof. 

This result illustrates one example where the Malthusian hypothesis, under we ex¬ 
amined the growth-fragmentation equation, fails, and where the stochastic model X 
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(whose mean measure could otherwise be expected to give a solution) does not remain 
locally finite. However, since uniqueness generally fails for the self-similar growth- 
fragmentation equation, it does not immediately imply that there is no global solution 
of (1.2). 

The procedure of creating the growth-fragmentation process y can be carried out 
under the general conditions on a, b and K given in the main body of the paper (this 
is done in Bertoin [3]) and the self-similar time-change (5.1) can also be applied in 
this general context for any a yf 0 (see [4, Corollary 2].) However, it remains an open 
problem to determine necessary and sufficient conditions for the process X to be locally 
finite at all times, and to decide when global solutions of (1.2) exist. 

6. Branching particle system and many-to-one formulas 

In this concluding section, we aim to clarify the connection between our work and 
the ‘spine’ or ‘tagged fragment’ approach to branching particle systems and fragmen¬ 
tation processes, with the hope of explaining the source of the solutions obtained in 
Theorem 3.1 and Corollaries 4.1 and 4.2. We refer the reader to the survey of Hardy 
and Harris [15] for results in the context of branching processes, Bertoin [2, §3.2.2] for 
the background on fragmentation processes, and Haas [13] for the use of the tagged 
fragment in solving the classical fragmentation equation. For the sake of simplicity, we 
focus on the case when the dislocation measure K is finite and the operator £ has the 
form (1.3). 

Let us assume that we can construct a system of branching particles in (0, oo), with the 
following dynamics: particles evolve independently of one other, each particle located 
at a; > 0 grows at rate cx°‘^^, and a particle located at a: > 0 is replaced by two particles 
located respectively at xy and a:(l — y) at rate x°‘K{(ly). Let Z(t) = {Zi{t))i>i denote 
the collection of particles in the system at time t > 0, starting at time 0 from a 
single particle located at 1. Informally, the verbal description of the dynamics of the 
particle system suggests that for every test function / G C“(0, oo), the functional 
F’(z) = for z = belongs to the domain of the infinitesimal generator 

G of the process {Z{t))t>o, and that 

= ^z?‘(czif'{z^) 4 - [ if{y I Zi) - f{z^))K{dy)). 

,>i V d[l/2,l) / 

Therefore, if we write Ht for the intensity measure of Z(t), i.e. 

(Mt, /) = E(£(Z(t)) = fife fiZ^))], 

^ i ■' 

then Kolmogorov’s forward equation entails that {pit)t>o solves the fragmentation 
equation (1.2). 

The analysis of the system (Z(t))t>o can be significantly simplified by identifying a 
spine among the particles, and formulating questions about the entire system in terms 
of just the spine particle via a many-to-one formula. This proceeds roughly as follows. 
Suppose that we can identify a function (t, z) Lp{t, z) > 0 such that the process 
M{t) = X]i>i Zi{t)) is a martingale. We introduce a new probability measure P 
by means of a martingale change of measure using M, simultaneously identifying one 
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of the particles to be the spine; specifically, we identify particle i of Z(t) as the spine 
with P-probability proportional to Zi{t)). 

At each time 0 < s < t, the spine particle has a unique ancestor, and we define the 
random variable W (s) to be the position of this ancestor. We now aim to identify the 
law of certain functionals of Z in terms of the law of W. In particular, if we define pt 
to be the law of W{t), we obtain 

(/r„/)=E(|]/(Z,(t))) =(p„/Mt,.)), 

which is known as a many-to-one formula. 

The spine method for solving (1.2) can be summarised as follows. We first use the 
dynamics of the branching particle system and the effect of the martingale change of 
measure to identify the process W. The one-dimensional distributions of W give the 
collection of measures {pt)t>o, and then the many-to-one formula gives us an explicit 
description of 

The method we have sketched can be made rigorous in the homogeneous case a = 0, 
even in the more general situation when the dislocation measure K is infinite and fulfills 
(1.5). More precisely, one can take (p{t,z) = exp{—K{uj)t)z‘^ for any w G domK, and 
then the process W(t) is the exponential of a Levy process with no positive jumps and 
Laplace exponent = k{oj -!-•)“ '^(‘^); this is the justification for Remark 3.2. We 
stress, however, that the general self-similar case a 0 is far less simple. In particular, 
it is not clear whether the branching particle system can indeed be constructed since, 
as noted in the previous section, explosion may occur. A fairly general class of growth- 
fragmentation processes was introduced recently in [4] by means of a Crump-Mode- 
Jagers process; however, although it is expected to be related to growth-fragmentation 
equations as described above, so far no many-to-one formula is known to make the 
connection rigorous. 
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